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Abstract 

Starting from Jefimenko's equations, we consider the multipole expansions of electric and mag- 
netic fields for a confined system of charges and currents. We analyze and comment on the calculus 
of radiated power, on the consistent use of approximation criteria, on the invariance of physical 
results when changing the point of reference, as well as on the use of electric and magnetic moments 
described by symmetric and trace free tensors. 
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I. INTRODUCTION 



Despite the successful and long history of the electromagnetic field theory, there are 
several topics open to new theoretical and pedagogical contributions. One of them concerns 
the formalism of the multipole expansion of fields, in general, and of the radiated one, in 
particular. Another issue is related to the importance of Jefimenko's equations in the study 
of such problems. Motivated by recent publications on the topic, we discuss some features 
of this type of problems, the present paper being a revisited and completed version of a 
previous electronic preprint^. 

The multipole expansion of the electromagnetic field in Cartesian coordinates is exposed 
in electrodynamics textbooks, as the well-known Refs.- anc£. Ordinarily, these expansions 
are calculated only in the first or second approximation, the higher-order terms being consid- 
ered too complicated. As Jackson writes in his textbook, the labor involved in manipulating 
terms in the expansion of the vector potential becomes increasingly prohibitive as the expan- 
sion is extended beyond the electric quadrupole terms (see Ref.-, pp 415-416). For this reason 
and due to the applicability only in the long- wavelength range, another treatment, based on 
the spherical tensors and on the solutions of Helmholtz equation is preferred. This alterna- 
tive has also a larger domain of applications. Actually, starting from the results obtained 
employing this calculation technique, the reader can verify what effort is involved when 
returning to the multipole Cartesian moments which offer a higher physical transparency 
(see Ref.-). A relatively recent textbook- and a paper-, the last related to the importance 
of Jefimenko's equations for expressing the electric and magnetic field when discussing the 
radiation theory, brought our attention on the necessity of some supplementary explana- 
tions. There are some prescriptions in the literature^,-,- for calculating higher-order terms 
of the multipole series based on a simple algebraic formalism of tensorial analysis. One of 
the aims of the present paper is to show how one can hide, as much as possible, the higher- 
order tensors behind some vectors, reducing the calculation technique to the formalism of 
an ordinary vectorial algebra or analysis. 

We start in section [III by shortly presenting the notation convention we use by giving a 
general formalism for handling multipolar expansions in Cartesian coordinates. In section 
Hill we derive the electric and magnetic fields multipole expansions without using the re- 
tarded potentials, while in section HVl we present some fundamental ideas of expressing the 
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multipole moments by symmetric trace free tensors. In section |V] we discuss the approxi- 
mation criteria stressing the significance of <i/A-criterion. In section |VI| we further discuss 
some features of the calculation for the radiated power, with an emphasis on the 4 — th 
order approximation in d/X < 1. Section [VIII presents the problem of the translational 
invariance of the fields expansions. Finally, in appendix |XJ we give the guidelines for the 
general tensorial calculus of the electric and magnetic moments and some general results in 
terms of symmetric trace free tensors. The last section is reserved for conclusions. 



II. GENERAL FORMALISM 



Let us write Maxwell's equations with a "system free" notation (a notation independent 
of the unit system) : 

Vxfi=^(j + £o f), VxE = -If, V-B = 0, V-E = ±p (1) 

a \ at J a at e 

where e , p , a are proportional factors depending on the system of units and satisfying the 
equation 

— = c 2 . (2) 

c is the vacuum light speed. Maxwell equations written in SI units are obtained from 
equations ([TJ) for a = 1 and the SI values of e , fi . For the Gauss system of units, 

a = C, Eo = 1/47T, jlQ = An. 

The relations between fields and potentials are written as 

1 B A 

B = VxA, £ = -V$ — . (3) 

Ordinarily, for the field calculation when the charge and the current distributions are 
known as functions of r and t, one firstly determines the retarded potentials: 

v v 
The fields E, B are then derived from equations (jBJ). The domain V includes the supports of 
p and J, and R = r — r' . The solutions (T4j) are indeed electromagnetic potentials verifying 
the Lorenz condition 

a at 



This equation can be directly verified using the continuity equation at the retarded time 
T = t- R/c: 

^ + [V'J(r',Ok T = (6) 

and the relation between the retarded value of a spatial derivative of a function and the 
same spatial derivative of the retarded function: 

r) = [d[ f(r', t% l=T + ~f{^ r). (7) 

The last equation should be well-known for each student from a class of electrodynamics 
since when writing the retarded potentials as particular solutions of the wave equation, it 
is necessary to verify the Lorenz condition. Only in this way one can be convinced that the 
retarded potentials are indeed electromagnetic potentials. This verification can be realized 
in a direct calculation, a good exercise for the student. 

An alternative approach for the field calculation is given by Jefimenko's equations (see 
Ref.-, section 6.5). The fields are obtained directly as retarded solutions of the following 
wave equations which are consequences of Maxwell's equations: 

A „ 1 d 2 B fin Ar ^ 1 d 2 E 1 /i &7 



These equations allow the retarded solutions 



v, 



V 



V V 

Using equation ([7]), the following expressions of the solutions (Q can be obtained: 
/io w J(r',r) 3 , 



>P( r 'i T ) j3 / f j(r',r) , 



p 

E(r,t) = — / V^-^dV--^ / " v " ,,J dV, 

V V 

where the dot signifies the time derivative. These last expressions can be obtained also 
introducing equations (00) in equations ([3]) and inverting the order of the derivatives and 
integrals. Considered as solutions of equations dHJ, the expressions (fTUI) are known as Je- 
fimenko's equations. We think it is important to show that sometimes it is convenient to 



perform a calculation directly on the fields represented by equations (fTO]) instead of calcu- 
lating firstly the potentials and after that, applying equations Q. In fact, in the present 
work we benefit from this circumstance. We owe our inspiration to Ref.-. 



III. MULTIPOLE EXPANSIONS OF FIELDS 

Let us consider the field B(r,t), given by equation ffTU]) . in a point from the exterior of 
the domain T> which includes the supports of p and J: 



\r — £| 



dV. (ll) 



V V 

Ei are the unit vectors of the Cartesian axes. Writing the Taylor series about £ = of the 
integrand, we can finally write 

B(r, t) = £- e i£ijk [ dV E ^ < - < *A ■ ■ ■ ^^f^ (12) 

where To = t — r/c. Here, we can apply a formula for the multiple derivative of a function 
of the type F(r )/r: 

U h--- U in /—I c n-l r l+l h-in Q t n-l ' ^ 10 ' 

1=0 

The coefficients are symmetric in ii, . . . , i n and can be expressed as 
[f] 

cfc.2. = • • -^-i^^h-i • • • (14) 

fc=0 

In the last equation, [f3] is the integer part of j3 and u = r/r. By A^...^} we understand the 
sum over all the permutations of the symbols i q giving distinct terms. For the objective of 
the present paper (n < 3), the coefficients D in equation ( fT4|) can be calculated directly from 
the successive derivative operations but, generally, one can establish reccurence relations. 
We can write now the expansion (112!) as 

B(r t) _ ^0 u V ('I)" V 1 C (n+1,0 /" d 3 / / / ^"^(tVq) 

n>0 ' £=0 ^ 

Next, we admit that in the previous equation one can invert the derivation with the integral 
operation on the domain T> and we have 

,. (_-\ \n I nn+l-Z 

B(r>t) = £b<^E^r- E^n^HT^s^^-*'^' (16) 

n>0 ' Z=0 



where it is introduced the n — th order magnetic moment fJJl^ with the Cartesian compo- 
nents: 

m ..,, U) J x' il ...x' in _ l J in (t)d 3 x'. (17) 

V 



From equation ( TTOI) . performing a similar calculation as in the magnetic field case, one 
obtains the multipole expansion of the electric field: 

1 ( 1 \n n+1 1 an+l-l 

Jt,[r,t) - --n— e i2_^ I 2^ c n+l-l r l+l U i<i...«n Q+n+l-i ^i-irA 7 ^ 
U n>0 ' i=0 



re; 



47ra 2 ^ n! ^ c n-; r «+i 

n>0 Z=0 

Pii...i„ are the Cartesian components of the electric moment of order n, P^: 

Pii...in(*) = / x' li ...x' ln p(r',t)d 3 x'. (19) 
v 

We point out that the expressions from equations (fTBjl and ( fT8l) are obtained without using 
commuting properties of the different operations except for the integration and the Taylor 
series expansion and, also, for the time derivation and the spatial integration. Of course, all 
the necessary convergence properties are supposed satisfied. In equation ffT6l) and f|T8|) one 
recognizes the multiple spatial derivative of 9Ri 1 ...i n k('r )/r and Pi 1 ...i n ( T o) / r as, for example, 

"+ 1 i Xn+l-l i 

E C^^^Tt^Ur.) = d 3 d h ... d ln [-^,., nfe (ro)] . (20) 

1=0 

The fields B and E can be expressed as 

B(r,t) = ^e i e&Y l ^d i a h ...a J^"- M 
Ana ' n\ r 

n>0 



i^vxVtffv-H^M (21) 

Ana 4^ n r 

n>0 



and 



Anen ^— ' n! r 47rcr ^— ' n! 

n>0 n>0 



(n+1) 



LvVt^V-II^M-^E^V-ll^— , (22) 



Anen t-r 1 n\ r Ana 2 n\ 

u n>0 n>0 
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where the following notation is introduced for the tensor contractions: 



Ail—in—mjl—jm-Bji—jm ■, n > m 



11— l| r 



A ■ R 



A . . R. . . . 



n = m 
n < m 



(23) 



We recognize, comparing these equations with equations (jl]), the multipole expansions of 
the potentials: 



A(r , f) = ^ E (z^y.,|2* ( " +1, M 



47TOJ 14 — ' 71! 
n>0 



47T£o " 
ri>0 



-l) n P (n) (Yn) 

-4-V n H- 1 ° ; . (24) 



Let us define 



*; C> = Cii...i„ y • • -4 n 4(^',i) dV = C»i...i„2Rii...*n*(*)j 



(25) 



with Cii...i„ symmetric in the n indices. Considering the following consequence of the conti- 
nuity equation 

, dp(r', r ) 



J k (r',T ) = V'-[x' k J(r',T )] + 



dt 



we can write 



a k {r,T ;(,n) = J x' k J{r',r ) ■ V{x' h ...i-Jd x' + C«...i„Pii..i„fc(7o)- 

Partial integrations and vanishing surface integrals are employed for obtaining the last result. 
Further, because of the symmetry of £, we get 

a k {r,T Q )(,n) - Cix...i n Ph...i n k(r ) = -n( h ... in / x' h . . . x' inl x' k J in (r', r ) dV 

v 

= -nCi!...i n £ki n g J . . . x- n _ 1 [r' x J(r',T )] q d 3 x' -na k (r,T ;(,n). 
v 

Let us introduce the standard magnetic moment of order n, , by its Cartesian 



components^ 
M,,.,,(/) 



n 



(n + l)a 



fx' il ...x' in _ 1 [r'xJ{r\t)\ in ^x'. 



(26) 



v 



We can write 



a k (r,T ;(,n) = 1 ...i n 9^ il ... infc (r ) = (h...i n -aejfei„,Mi 1 ...i B . 19 (r ) + — — P h ...i n k(r ) .(27) 



n+ 1 



This result inserted in equation ( |T6l) gives 



■1)' 



n-l 



n>l 



n+1 

E 

Z=0 



1 



_£,(n+l,0 



^n+1- 



/'() 



(1+1 



Ana 



,n+l— J r i+l J»i...in gpi+l-l 
C (n+l,l) d n+l ~ l p 



M il ... in _ 1 q(r ) 



(n + 1)! ^ cn+i-V^ 1 J"* 1 -*" n - in 

n>0 v y Z=0 



fc(ro). (28) 



With the same calculation technique, one obtains the multipole expansion of the electric 
field: 



E(r,t) 



1_ (_!)n ^ i 



£,(n+l,Z) ^ 



47T£o *^ n\ ^ C n+l-l r l+l Ul-in Q t n+l-l ^ 
u n>0 Z=0 



(to) 



/'o 



£»0M) 



& 



n—l 



(_l)n » ! 



Ana 2 ' ^ in + 1)! ^ C «-M +1 <9t n 

n>0 v ' 1=0 



Pii...i„i(To) 



/'() 



E 



■l) n 



4ira ' — ' n\ 

n>l 



E 



Ct 1 } I^M. 



'ii n q n -l l+l w il...i„ Q+n-l "U-Jn-KJ 

«=0 



(ro). 



(29) 



Processing the expressions from equations (128]) and fl29|) by the same procedure applied for 
obtaining equations (12~TT) and (|22|) . we obtain 



E 

n>l 



n! 



a ^ (n + 1)! 

n>0 v ' 



ly p (n+1) 



(30) 



and 



B(r,t) 



4vre 
a \ 

+ 



n>l 



^E 

n! 



-1)" „ M P W (ro) , 1 



n>0 
n-l 



??■! 



4)n 



P'" +I, (ro) 



(n + 1)! 

n>0 V ' 



V X V 



(31) 



From equations (130]) and (l3Tj) one can immediately identify the multipole expansions of 
the potentials A and $. 

The expansions from equations ( 1301) and (13T1) can be easily processed by a vectorial 
calculus introducing the following vectors: 



(32) 



8 



and 



r r 
Consequently, the expressions for the magnetic and electric field become 



B(r, ( ) = £vx 



and 
E(r,t) 



£ 

,n>l 



m-1 



1 



If 



111 



VxM(r,i;V,n) + -V / j — P(r, t; V, n + 1) 

a ^-^ (n + 1 ! 



(33) 



(34) 



n>0 



47T60 

a 



( 1 \n i / i \n 

v £ { —r- ( v • p ( r > *; v > ")) + -2 E A-ttt 7 ^ *; v ' ra + 

' n! c z *-^(n + l)i 

n>0 n>0 v 



« ("I)"" 1 „ „ x 

c * — ' n 

n>l 



(35) 



IV. EXPRESSING THE FIELDS BY SYMMETRIC TRACE-FREE TENSORS 



The general terms from the field expansions, given by equations (1281) and (1291) or (1301) 
and (13TI) . are apparently too complicated and they imply cumbersome calculations when the 
fields E and B are introduced in some physical interesting expressions as, for example, the 
radiation intensity. A considerable simplification of these calculations is obtained introduc- 
ing the irreducible electric and magnetic moments defined as symmetric trace- free ( "STF" ) 
Cartesian tensors. Let us consider a n — th order tensor and the corresponding projec- 
tions <S(T^) and T(T^ n ') on the subspaces of symmetric and STF tensors, respectively. 
For the electric moment P^ n \ we have a symmetric tensor and one has only to establish 
their trace free projection. Let us consider the simplest case of the quadrupolar moment 
p( 2 ). Writing the components as 



lj — < Hj 



T\n + AS. 



there is a unique value of the parameter A such that = T(P^). For A = P qq /S, 
1 _ „ f / 1 



n = p 

1 'i? r u 



P S- 

r qq U tJ 



(xiXj - -r 2 5ij) pd x. 



(36) 



v 



The STF projection of P*- 3 -* can be calculated searching the first order tensor A 1 - 1 ^ such that 
pj(3) _ X(P ( - 3 - ) ) is given by the components 



njj'fc Pijk ^{ij^k}- 



(37) 



From the condition of vanishing traces of the tensor n^ 3 \ one easily obtains: 



Ai = \P mi = \j r 2 x iP d 3 x. (38) 
v 

Concerning the magnetic quadrupole moment M^, we have a simple procedure for obtaining 
the STF projection. Let us write the identity 

where the first bracket represents the symmetric part of this tensor, and the second, the 
antisymmetric one. The symmetric part is, for this case (n = 2), a STF tensor = 
T(M {2) ). Therefore, 

My = Vij + - e ijk Njfe, (39) 



where 



N fc = e m My = A y [ r x (r x J)] fe d 3 x = J- y [(r ■ J)x k - r 2 J k ] d 3 x. (40) 
We consider now the effect of the substitution 

p(2)_> n (2). p..^ P .._ A5 .. (41 ) 

in the multipole expansion of the field. From equation (|3~2"1) we obtain the transformation of 
V(r,t; V, 2) by the substitution (gTJ: 

P(r, f; V, 2) P(2 ^ <2) P(r , t; V, 2) - V (42) 

r 



It is easy to see from equation ( 1341) that £? is invariant to the substitution ( 1421) . In the 
expression (|35p of the electric field E, the combination of the electric quadripolar terms 
from the first two sums gets 

A(t„) 1 Afro)' 



87T£ 



A 



c 2 r 







since, for r ^ 0, A(ro)/V is a solution of the homogeneous wave equation. Therefore, the 
field E is also invariant to the substitution 042p . 

The circumstances change dramatically for n > 3 in the cases of electric moments and for 
n > 2 in the magnetic ones. As we will see in the following, the electromagnetic field is not 
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invariant under the substitutions — ► n < - 3 - ) and M*- 2 - 1 — > r 1 - 2 ^. The general case, i.e. for 
rt > 3 for electric moments and n > 2 for the magnetic ones, is developed in earlier works 
(see Ref.- and the literature cited in). We insert in the appendix some general results from 
these works, without demonstrations. 
Let us consider the substitution 

p(3) ^ n (3) : P ijk P ijk - A {1 5 jkh (43) 

in the fields expressions. In this case, 

1 A(r ) 



P(r,t;V,3) - P(r,t;V,3)-2V 



r 



c 2 r 



P(r.,t;V,3)-2Vx(Vx^)-i^M, (44) 



where A = e^Aj. Performing the substitution (1441) in equation (134|) . we obtain the following 
transformation of B: 

B(r,t)-.B(r ) t)--^Vx^). (45) 
247racr r 

The same substitution in equation (I35I) gives 

J5(r, f) - J5(r, f) - V [V ■ AM] + (46) 



Both, the modifications of the magnetic and electric field given in equations ( 145]) and ( 1461 ). 
respectively, can be compensated by the following transformation of the electric dipolar 
moment: 

P-P+^A- (47) 
Let be the substitution 

M (2) ^r^: My-^My-^eyfcNfc, (48) 
corresponding to equation (1391) . Then, 

M(r, *; V, 2) - M(r, *; V, 2) - \ V x 

2 r 

where N = e^Nj, and 

V x M(r, t; V, 2) - V x M(r, i; V, 2) + \ V [V • ^] - ^ (49) 

11 



Performing the substitution (149]) in equations (1341) and (|35|) . we obtain the transformations: 

B(r,t)^B(r,t) + ^Vx^M. (50) 
Lb7rc z r 

and 

fi(r, t) E(r, t) + -BL V( V ■ * - N < T °> . (51) 

lo7ra r \bnac z r 

The changes of B and B described by equations (150]) and floTT) are compensated by the 

transformation of electric dipolar moment, too: 

P-P-^N (52) 

Finally, substituting in the field expressions the momenta P*- 3 - 1 and by their STF- 
projections \~\^ and V^ 2 \ respectively, the corresponding changes of the fields can be com- 
pensated by the transformation 

p-P = P-^(f *-±a) =p- it, (53) 

where by the vector t we understand 



t = — N 

4 6 



l A = ToJ ^ r - J ^ r ~ 2r2j )^ x - ( 54 ) 



v 



The last expression is obtained from equations (138]) and ([40]) applying the continuity equation 
together with an operation of partial integration. This is the so-called electric toroidal dipolar 
moment, its presence in equation ( |53]) leading to physical effects similar to the ones produced 
by an electric dipole moment. This fact was observed firstly in Ref.— (see also Refs.— ,—,—)). 

V. APPROXIMATION CRITERIA 

Generally, one can operate with the spectral decomposition of the charge and current 
distributions: 

oo oo 

p(r,t)= j A.We-^dw, J(r,t)= J J w (r) e"^ du. (55) 

— oo — oo 

Correspondingly, one can introduce the spectral decomposition of the multipole moments: 

oo oo 

p(»)(t) = y pW e - w do;, M (n) (t)= J M^e-^dw. (56) 

— oo — oo 
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The Fourier components are given by 

= /"<••• < Pc(r') dV, M Hil ... in = / x' h . . . x' ln _ x (r' x J^{r')). dV. (57) 



v v 
Let us denote by d the linear dimension of the domain T>. For estimating the order of 

magnitude of the different terms from the multipolar expansions of the fields, we notice the 

following relations: 

P (n) ~ d n , m {n) ~ d n+1 , M (n) ~ d n+1 . (58) 

Considering the Fourier components in equations (1161) and (ITS]) , we can estimate the con- 
tributions of the parameters d, A = 2-kc/uj and r to the orders of magnitude of the different 
terms from the series expansions. For example, the estimation of the order of magnitude in 
equation (TT6I) leads to the relations: 

n+1 -, pm+l—l 

c( n ) (, ,\ _ \^ 1 Mn+1,1) O / -\ut+\ujr/c\ 

°ji 1 ...i n k\ UJ ) Z-^i c n+l-l r l+l ^jil-in ftm+l-l \ MV {u)n...in-ikV ) 

1=0 

l^ U n+l-l d n+l x ^ d n +1 

Y ' 4 qH+1—1 y J ^^+1^' fl 

1=0 1=0 

that is 

n+1 



s (n) ■ ,M ~ - 

, - y jii...i„fcV Lv ; 



A /A xn+1 
1 + - + ... - 
r V r 



(59) 



From equation (1591) we can conclude that for the convergence of the expansion series it is 
necessary to have the inequality 

A > d (60) 

satisfied. Consequently, for the spectral decomposition (155|) . the extension must be limited 
to the wave lengths verifying the last inequality. We can conclude that there are three 
distinct regions around an isolate source with corresponding approximations of fields: the 
near (static) region d < r << A, the intermediate (induction) region r pa A and the far 
(radiation) region r >> A, Ref.-. We consider here only the last region. 



VI. THE RADIATION FIELD 



For r >> A, the dominant term in the bracket from equation (1591) is represented by 
the unity. The corresponding fields E and B are proportional to 1/r and represent the 
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parts contributing, as it is well known, to the energy and linear momentum radiated at 
large distances. The next term from the same bracket corresponds to the part of the field 
proportional to 1/r 2 which contributes to the radiated angular momentum. In the present 
paper, as in most textbooks of electrodynamics, we consider only the radiated energy and 
therefore, we calculate the part of the field proportional to 1/r. However, we point out 
that, as basic variables, the fields E and B must be defined such that all observables of the 
radiated field can be calculated. The angular momentum is one of these observables (for a 
complete calculation of the radiation field, see Refs.-,- and literature cited in). 

Let us calculate the first dominant terms, proportional to 1/r, resulting from equations 
( ITS"]) and (ITS]) . For this purpose, we retain only the terms corresponding to the coefficients 

c£Z = (-l) n v il ...v in . (61) 
The expression from equation (Tl6|) becomes 

mo 1 dn+1 

n>0 

or 



B -<^ = ; " X E |sT (""I |3*<» +1 >(t„)) (62) 



B rad (r, t) = --^ z J2 ( " x ^TT m ( r o; "> n + X ) ) > ( 63 ) 



4ira r 

n>0 



where we introduced the vector 

m(t;£,n) = ^ Tl ~ 1 | |SDT <Tl) (*). (64) 
Similarly, equation ffl~8j) gets 

n>l 

fi 1 ^ 1 9 



I y _L ( v n\ ^(n+i) ( <\ (65) 

Ana 2 r ^ n!c n <9t n+1 V 11 v U V v 7 



n>0 



since the source is an electric neutral system (P^ = 0), or 

1 1 f) Tl + 1 

n>l 

/i lr- 1 d n+1 



Ana 2 r ^— ' n!c n <9t n+1 

n>0 
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with 



7z(t;tn)=e^\\P {n \t). 
Based on equation (1271) . we can write the following relation for the vector m: 



m(r ; u, n + 1) = -av x h(t ; is, n) H —t^(t ; u,n+l), 

n + 1 



where 



(67) 



(68) 



(69) 



Substituting equation (|68|) in equation (163|) . one gets 

B rad (r,t) = ^ - <j . Li L +1 [" x x M(^o; v,n))} 



An r I ^ n!c n+1 <9t r 

. n>l 



1 1 d 

- « E (n+ l) !c n + l^ +2 [» >< - + 1)] 

The change n — > n + 1 in the summation index of the second sum allows us to write: 



(70) 



B rad (r,t) 



/io 1 1 d n+1 

V X > 

An r ^ n\ c n+1 dt n+1 . 

n>l 



v x fi(r ; v, n) 7r(r ; v, n) 



a 



(71) 



Applying the same procedure to equation (166|) . we obtain 
E rad (r,t) 



1 1 V- 1 d n+1 \, , u a . 



n] 



l r 1 «9 n+2 , 1N 1 
c ^ (n + 1 !c™ +1 dt n+2 v ' ' ; f 

n>'0 V 1 ) 



(72) 



The same change of the summation index n — > n + 1 in the last sum gives: 



E rad (r,t) 



1 1 



-v x 



1 <9 



■n+l 



4vt5 r ^ n\c n+l dt n+1 L 

n>l 



i/ x 7r(r ; u, n) H /x(r ; v, n 



(73) 



Comparing equations (17T1) and (1731 . we notice the relation: 

E rad [r,t) < — > B rad [r,t). 

As seen from equations (1711 and ( 1731) . the fields -E ra <i and B rad are purely transverse fields 
satisfying the properties: 



v ■ E rad — 0, f • -B ra rf — 0, E rad — — B rad x v, Eo\E rad \ — — \B rad \ 

a no 



(74) 
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i.e. the characteristic wave properties. 

Considering the wave region, the essential parameter characterizing the multipole expan- 
sion is d/X. Coming back to the monochromatic case studied in Section |VJ and introducing 
the notation C = f << 1, we can conclude that 

dt k \ C k , k<n dt k | C fc+i 5 k < n 

In the field multipolar expansions in terms of primitive moments P^ and as, for 

example, in equations ( ITUl) and ( 1721) . the series terms appear in an increasing order in the 

~ (n) 

parameter (. If the multipolar series are expressed in terms of the STF-tensors M and 

~(n) 

P as, for example, in equations (1701) and (I72|) . some caution is necessary. Replacing the 
pair M^ -1 ), P^ by the corresponding STF-projections, one induces in the moments of 
the inferior order compensating terms of orders ( k , with k < N. If we are interested in an 
approximation which includes the multipole moments up to a given rank n < N, we must 

~(m) ~ (m) , 

discard the compensating terms from P , M , m < n of orders ( for k > n. 

The task of writing finite multipolar sums in a well defined approximation becomes del- 
icate when one considers products of field quantities as, for example, when defining the 
radiation intensities. The simplest, and most important, is the case of the radiated power. 
The Poynting vector can be written in terms of E ra d and Brad with the help of equations 
(G3D as: 

S = —(E x B) = e \E rad \ 2 cv + 0(1) = -\B rad \ 2 cv + 0(1). (76) 
Ho r A no H 

B and E are considered real vectors. The total radiated power may be written as the limit 

of a surface integral on a sphere centered in O, of radius r, for r — > oo. Let us express the 

energy current in the radiation approximation corresponding to the sphere of radius r: 

N{S,E r ;t)= I r 2 u ■ S(r,t)dn(u) = — I r 2 \B{r, t)\ 2 dfi(i/) + 0{-). (77) 

The integrand from the last expression represents the angular distribution of the radiation. 
The quantity N(S, S r ; t) dt represents the energy which crosses the sphere S r in the time 
interval (t, t + dt) and is determined by the values of the multipole moments of the source 
in the interval (t — r/c, t + dt — r/c). For large, but finite r, r >> A, this is an approximate 
expression obtained by neglecting the terms of order at least 1/r. The expression from 
equation (I77|) can be employed for drawing conclusions on the electric charge distribution 
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at the retarded time from observations, on the angular distribution or on the total radiated 
power. Since dt = dro, we can say that N(S, S r ; t) dt is the part of the energy emitted by 
the source in the given time interval which contributes to the energy intensity correspond- 
ing to S r at the time t + r/c. If in equation (1771) we put t instead r , then the quantity 
lim^oo N(S, £ r ; t) dt represents that part of the energy emitted by the source which con- 
tributes to the radiated energy or, shortly, radiated by the source. The situation changes 
when the support of the source depends on time, in particular, for the radiation of a moving 
poit-like source (see Ref.-,&73). 

For calculating the angular distribution of the energy radiated or the total radiated 
power, one deals with the square of B rac i or E ra d. Let us consider firstly that we operate 
with the series (170j) or (172"j) expressed in terms of primitive electric and magnetic moments 
. For illustrating the kind of problems which appear in the approximation process, it is 
sufficient to consider the well known problem of the calculation of the electric and magnetic 
dipole and electric quadrupole field contributions to the radiation (see Refs.-,-). From the 
very beginning, we point out that here, this calculation is associated with the problem of 
the radiation of a complex system characterized by an arbitrary number of electric and 
magnetic multipoles. Once considered the contribution of the electric quadrupole field, from 
the expression of \B ra d\ 2 , with B considered real, the square of the electric quadrupolar 
moment is part of this contribution. This is a term of the order ( 4 . There is a similar 
conclusion concerning the contribution of the magnetic dipolar moment m. It is easy to see 
that in a consistent approximation procedure, together with the squares of P^ 2 \ m and p, we 
must consider also the products of the electric dipolar moment p with the electric octopolar 
moment and the magnetic quadrupolar moment M^ 2 K They also give contributions of 
the order ( A . Let us estimate these contributions in terms of primitive moments. In equation 
(|76|) . to the contribution of (Bf^ d + B^ d ) 2 , we have to add the fourth-order contribution of 
the product B^] d ■ B^ d . From equation ( !70l) we can write, for n = 0, 1, 2: 




IMi 1 
47rac r 



v x 7t(t ; v, 1), 



(78) 




Anc 3 2r 



Ho 1 

4 7rc 2 r 




v x \y x jj,(r ; v, 1 



v x [i/x fi (t ; v, 2 



)] 



1 



1 



ux 7T (r ;i/,2) 



i/x 7T (r ; ^,3) 




(79) 



(80) 
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The angular distribution and the total power of the radiation are given by 

X {u) = -r 2 {B rad )\ 1 = ^( r "(B rad ) 2 }, (81) 

where 

(82) 



(/(«')> = ^ y /M^M- 

Considering, for example, the angular distribution X(i/), expanded up to the order ( 4 , we 
have 



X(u) = f r 2 + B« ) 2 + 2B^ ■ S3 + 0((% 



.(0) n (2) 



/'() 



(83) 



In the following, for simplifying the notation, we understand by 7r(n) and /x(n) the corre- 
sponding vectors with the arguments To and v. A simple algebraic calculation gets: 



±-2/1,(0) , d(1)\ 2 = 1 

/X 1 ^ ^ " (47T) 2 £ C 3 



+ • ("X (2)) + (2)) - - C (U ■ TT(l)) („• 7T (2)) 



(84) 



(4tt)% c 5 



a ji (2)- (i/x *(!)) + -(*(!)■ £(3)) 



--(i/.7r(l))(i/.5?(3))j. (85) 
One can determine the averaged quantities from the last two equations using formula^: 



0. 



n = 2k + l, 
n = 2k, A; = 0,1, 



(2fc+l)H ^{*i*2 • • • Sin-iin} 

The averaged expressions from equations (184"]) and (185]) are given by 

1 



± -2/ / R (0) , d(1) ^| 2 \ _ L 

{ [ rad + rad) ' ~ 47T£ 0C 3 V3 " ' 3C 2 



2 „ 2 2a 2 .. , 



20c 



2 ■ ij • ij 



60 c 



2 *» 



and 



2c 



\^rad ^ rad) 



1 



/i 47re c 5 V 3 45 



(86) 



(87) 



(88) 



In equations (IH7]) and ([881 some standard notation for magnetic and electric moments is 
used. A simple calculation shows that 



— r 2 (B (0) -B {2) )- 1 4 ft t 

/i r * l '~ 4VT5 0C 3 3C 2P ' 



(89) 
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where t is defined by equation (|54|) . The total radiated power is then given by 

1 (2 „ 2a 2 „ 1 - - 1 - 2 4 -\ , . 

Z = ? - p 2 + — m 2 + P i? P i? - P.. p ■ t . 90 

4ne c 3 V3 3c 2 20c 2 J J 60 c 2 3c y V ' 

Considering in equation ( 1901) the multipolar moments at the time t, we have the description 
of the source emission. This expression differs from the result given in Ref.- firstly by the 
fact that, whereas here the electric 4-polar moment is the primitive one, in this reference it 
is represented by the corresponding STF tensor and, consequently, the term containing the 
trace of the electric 4-polar moment is absent. The second difference consists in the presence 
in equation fl90l) of the toroidal moment t. 

For expressing X by STF- tensors, we apply the invariance of the field to the substitutions 
of all multipole primitive tensors P and M by corresponding STF-tensors P, M . These last 
tensors are, in general, different from the correspondent STF projections as it was shown in 
Section [IVl Obviously, for the approximation corresponding to equation (190]) . it is sufficient 
to consider the STF-projections P = T(P (3) ) and M = T(M (2) ) and the induced 
transformations — ► ^ for k — 1, 2, and M 1 -^ — ► ^ for k — 1 (see Section HVl 
equations ( I43|) -(l54l)). Performing these substitutions in equations ( 1H71) and ( |88fl . we obtain 
for X: 

The term t 2 /c 4 must be eliminated since it is of order £ 6 , such that the final expression of 
X in the required approximation is given by 

1 (2 .. 2 2a 2 .. 2 1 4 .. -\ 

This is the correct expression of the radiated power in the ^-approximation and it is obvi- 
ously not represented by independent contributions of p, m, V\^ 2 \ In Ref.-, the contribution 
of the multipole moment P*- 2 -* is calculated as an independent contribution. The correspond- 
ing result is a correct one if the radiation source is an elementary system characterized only 
by his electric 4-polar moment. The reader must be warned that in case one considers a 
complex system characterized by higher multipole moments, the correct evaluation of the 
contributions of the moments p, m and P^ cannot be identified with the sum of the inde- 
pendent contributions of these multipole moments. Particularly, this is the reason why in 
Ref.- (& 71) one obtains an incomplete result when calculating the electric and magnetic 
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dipolar and electric quadripolar radiation considered as a result of the multipole expansion 
of a complex system. 

A simpler procedure for obtaining the expression from equation (I92p is achieved if 
in equations f[M|) and ()85[) one substitutes ir(t;i/,n) — > 7r(t;u,n) = i/ n-1 ||P (t) and 
n(t;u,n) — > fj,(t;u,n) = i/^HM (t). This procedure has the advantage of a simpler 
calculation of the z/-averaged quantities for higher order terms of X. In appendix |A] we give 
some general results for the radiated energy and linear and angular momenta (for detailed 
calculations see Ref.- and the issues cited there). 



VII. TRANSLATION AL INVARIANCE 



Let us consider a new point of reference O' specified in the Cartesian reference system 
{O, (ej)j =1 _;_3} by the vector a = aj e^. The new reference system in the affine space is defined 
as {O', (ej)j^_3}. Here, both systems are associated with the same basis in the vectorial space. 
For avoiding a possible confusion produced by the notation used, we denote by P the point 
where the field is calculated and by P' the current integration point. Correspondingly, 
instead of B(r,t) and Jk(r',r) in equation fTTT]) . the notation B(P,t) and Jk(P',r) will be 
used. Generally, the multipole moments depend on the point of reference chosen. As it is 
stressed in Ref.—, the "origin dependence" of multipole moments forms only a part of the 
effect of changing the point of reference. Considering consistently all the the changes in 
the calculation of an observable of the physical system, one can verify the invariance of this 
observable. We can write equation fTTT]) in the new system of reference as 



^"'^/ K-iT' ) dV ' (93) 

& \ ' S ' U=r" 

where r" = x" is the vector associated with the pair O' P' . Writing the Taylor series about 
£ = and, after that, performing the substitution £ = r", one gets 



td'i n jA ^° i~ 1)™ f a a a a a Jk(P',t ' ) 3 

B{P,t) = - — eiEijk^ — / x h ...x in djd h ...d ln d 6 x 

n>0 



p 

T,(pi + - l r ~ Q K 

3 / 



^0 ( _1 ) n l( i \ (i \a a a J k(P',t - JL ~^ j 

eiS ljk j— / ( x h - ah)--- {%i n ~ <knM d h ■ ■ ■ d in , , — d x 



Ana n\ J \r — a 

n>0 p 
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or 



47ra ^— ' n! ^— ' ra!m — m 

n>0 m=0 v 

/a a a / / Jk(P ) ^ — c ' ) ,3 / 
Oi Oi, ...Oi X 4 , . . . X, : r- 5 QI . 

V 

Performing the Taylor expansion of the fraction Jk(P',t — \r — a\/c)/\r — a\ as a function 
of a about a = 0, it results: 

-IV 



/ 1 \ 772 



47ra / -— ' ^—*m\(n — m)\ 1 m / -— ' /! 

n>0 m=0 v y Z>0 



• • • a jl 



xdj d n ... d in d n . . . d n ^ J x' m+1 . . . <J fe (r', t - r -)j dV. (94) 

This equation can be written, in a more compact form, as: 
B'(r,t) 

- ^ be B V V ( " 1)W ~ m (a vrV^'fa VV /V™ U ^222M 



Ana 1 lJ * 3 ^^m\(n-m)\ y ' ^ l\ 

n>0 m=0 v * l>0 



k 

= JpLe i e ijk d j 1 £ib Tj ( a ' V) m e~ a v f V"-| | g^^Z°A . (95) 

n>0 m=0 v ' \ / k 

Obviously, considering the infinite series corresponding to the multipole expansion of the 
magnetic field, we have no dependence of the choice of the reference point since the Taylor 
series gives the same values of a function no matter what is the point about the expansion 
is taken. Indeed, considering the last expression of B', for the sums after n and m, we have: 

oo n oo oo oo oo 

B '( r > t) = ^2^2fnm = ^2^2^(n- m)f nm = ^2^2 fnm, 
n=0 m=0 n=0 m=0 m=0 n=m 

i.e. the possibility of inverting the order of the summation after n and m in the case of 
infinite series. Using this property in equation fl95l) and introducing a new summation index 
q = n — m, we obtain 

Recognizing in the last expression the series of exp(a ■ V) and the expansion fl2TT) of B(r, t), 
one can conclude that B' = B, i.e. the translation invariance of the magnetic field multipole 
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expansion. For practical calculations, we have to take into account the expression in equation 
(195]) as a double series expansion upon the parameters a and d. Since a < d, the parameters 
a and d are associated with the subunit and adimensional parameter ( = d/X ~ a/ A. 
Obviously, the translation invariance cannot be verified by each term of the field expansion 
but, for a given approximation defined by a maximum power of the parameter (, it is verified 
for a well defined partial sum from the series. For illustrating how this invariance property 
is working, we write equation (1951) as 

B'(r,t)=j:±^^e-^B^(r,t). 

n>0 m=0 

Let us consider the first three terms from the field series: 

B '(0) = e -o.V B (0) = h _ ( a . V ) + ^(a • V) 2 - -(a ■ V) 3 + . . . 1 B (0) 

2 6 

B'W = [(a ■ V) - (a ■ V) 2 + i(a • V) 3 + . . . ] B (0) + [l - (a • V) + ^(a • V) 2 + . . 

£?'( 2 ) = [i(a- V) 2 -i(a- V) 3 + ...]B(°)+ [(a- V) - (a ■ V) 2 + ...]£?« 
+ [!_ ( a . V) + ...]# {2) . 

Since (a- V) is associated to the parameter ( and B^ and B < - 2 ' ) are of the orders (, ( 2 
and £ 3 , respectively, we can express the partial sums up to the order ( 4 inclusively: 

B '(0) + B '(l) = B (0) +B (l) + [_I( a . V) 2 + I (a . V )3] B (0) 

2 3 

+ [_(a-V) + i(a-V) 2 ] J BW + 0(C 5 ), (96) 

B '(0) + B 'W + B '(2) = B (0) + B (l) + B (2) 

- ha ■ V) 3 B (0) - ha ■ V) 2 £ (1) - (a • V)B (2) + C(C 5 ). (97) 
6 2 

From equations (1961) and (1971) it is seen that the sum B^ + B^ is translational invariant 
up to the order ( 2 and the sum B^ + B^ + B^ 2 \ up to the order ( 3 . An interesting 
circumstance appears in case one analyzes the second powers of these sums as, for example, 
when one calculates the radiation intensities. Let us consider the squared sum from equation 
©: 

( B '(0) + B '(D) 2 = ( B (0) + B (D)2 + (9g) 
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Terms of the orders ( 3 and £ 4 are present in (B^ + such that we have to ensure the 

translational invariance up the order four in the parameter (. Therefore, the contribution 
of B ^ must be added up since this term introduces the product B^ ■ B^ ~ £ 4 : 

( B '(o) + B 'd) + B 'my = ^(o) + B (i)j2 + 2jB ( ) . B(2 ) + (gg) 

In terms of the multipole moments 971, we conclude that, if one includes the square of 
the magnetic moment then, for ensuring the translation invariance, the inclusion of a 

partial contribution of Wt^ is necessary, too. This partial contribution is expressed in terms 
of the contraction of the tensors 971^ and SPT*- -'. Correspondingly, including the squares of 
the electric 4-polar moment and of the magnetic dipolar moment m we have to include 
also partial contributions of the moments M*- 2 -* and P^ 3 -* represented by contractions of the 
corresponding tensors with the electric dipolar moment p. This property explains why in 
Ref.— , introducing in the expression of the radiated power only the sum of the moments p, 
m and P^ squared, finally, for establishing a translation invariant result one performs a 
correction by a term represented by contractions of and P^ with p. In such a way one 
introduces in fact a contribution of the dipole toroidal moment t. 

The procedure can be continued for higher orders. In the following step, for example, the 
contribution of the square of B^ + B^ + B^ up to the order ( 6 is obtained as 

( B '(o) + B '(D + B'my = ^(o) + B (i) + B my + 2 ( B (o) . B o) + B d) . B (3)) (100) 

Besides the contributions of the multipole moments p, m, P^ 2 \ M^ 2 \ P^ 3 \ we have to include 
also some partial contributions from P^ 4 -* and represented by the last two terms from 
equation fllOOp . 

Combining the above observations with the results of the previous section, we can con- 
clude here that all the invariance properties are ensured if a consistent approximation pro- 
cedure is used. For this reason, we suspect some inconsistencies in considering the approx- 
imation criteria when, performing multipole expansion, one concludes that a correction is 
necessary for ensuring the translational invariance (see Refs.— ,—,-). 

VIII. CONCLUSION 

Along this paper, we mainly exposed the fundamental grounds for the calculus procedure 
of the multipole expansion of the electromagnetic field. 
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Firstly, the direct approach in the expansion of the quantities E and B from Jefimenko's 
equations revealed that this method is neither more difficult, nor easier than the one based 
on the expansion of the potentials A and $. An advantage of the current approach can 
be the elimination, at least partly, of the difficulties related to the inversion of different 
operations, mostly of derivation and integration. Another advantage is offered by the direct 
examination of the invariance of the fields E and B to the transformations required by the 
manipulation of the electric and magnetic moment tensors (sections II III and IIVj) . In this 
case, one does not need to employ gauge invariance. 

Secondly, in section |Vl we specify the approximation criteria according to which one 
has to consider three approximation intervals for the field associated to a localized source, 
as mentioned in Refs.- and — . Then, in the next two sections, for radiation calculus, we 
consistently apply the criterion d/X. The consequences are the ones mainly known from 
RefA 

We hope that the discussion from section IVIII regarding the invariance of the physical 
results when changing the point of reference, although without an extension to the general 
multipole terms, is convincing enough. We concluded that the translational invariance is 
automatically ensured if one respects the approximation criteria when employing multipole 
expansions. 

Since in the main part of the current work we wanted to avoid the generalization to 
superior orders of the multipole series, in the Appendix we summarized for the interested 
reader, the main known results in the literature. 

We believe that in an electrodynamics lecture, there should exist some remarks on the 
contributions of the toroidal moments, at least for the quadrupole electric moment. In this 
approximation, as one can see from sections IIVI and IVIl the highlight of the contributions is 
not involving a considerable computational effort. 

APPENDIX A: SOME BASIC FORMULAE FOR GENERALIZING THE PRO- 
CEDURE OF REDUCING THE ELECTRIC AND MAGNETIC MOMENTUM 
TENSORS 

In this section we list the required formulae for the generalized tensor reducing procedure. 
The equations are given without the related proofs since they can be found in the literature 
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(see Ref.-). 

Let be a symmetric tensor of rank n. Its STF- projection results from the following 
formula: 

T(SW). = S,,., n - 5 {!l!2 A(SW) t3 ..,„ } . (Al) 

The operator A is defined on the account of a formula for the STF-projection of a symmetric 
tensor given in Ref.- (the book^ is cited as the origin of this formula): 

A(s(n)) _ [n ^ ] (-ir[2n-l-2(m + l)]!! 

^ '<i-in-2 (m+ l)(2n- 1)!! 

X 5{i li2 . . . ^2m-li2m S i 2 ^ + t..i n _ 2 }- ( A2 ) 

The proof can be found in Ref.—. The notation S^ n ' p ^ indicates p pairs of contracted indices. 
In the following, for simplifying the notation, all arguments of the operator A should be 
considered symmetric tensors, i.e. A(T^) = A(«S(T'™' ) )) for any tensor ~P n \ The same 
applies to the operator T: T(T^) = T"(«S(T^)), «S being the symmetrization operator. 

In the symmetrization process we have to calculate the symmetric projections of some 
tensors of the magnetic moment type: they are symmetric in the first n — 1 indices and 
the contraction of i n with any index i q , q = 1 . . . n — 1 gives a null result. For the symmetric 
projection of such a tensor, we introduce the formula: 

1 n— 1 

( 5 ( L(n) ))„... = U-in - -E W^,„_ 1? (L (n) )> (A3) 

A=l 

where A/"-^ $ is the component with the index suppressed. The operator A/" defines 
a correspondence between and a tensor of rank [n — 1) of the same type from the 
symmetry point of view. Particularly, 

Af2k, M V)\ = (-!) fc rc /" (r2) /c r n-2/c x Jd 3 

v 7 (n + l)a y v ; 

^2/c+l ^(n)) = iZ^L J ( r 2)fc r «-2fc-l x (r X J) d 3 X, fc = 0, 1, 2, . . . , (A4) 

V 

where a n x 6 is the tensor defined by the components {a n x b)^ . = . . . aj n _ 1 (a x i>) in . 
Let us consider the process of STF reducing for the multipole tensors starting with the rank 
N + 1 in the electric case, and in the magnetic one. A general formula for the resulting 
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~(n) ~(n) 

tensors P , for n = 1, . . . , N and M , for n = 1, . . . , TV — 1 is given below. It includes 
all compensating terms obtained in this process 



p (n) = r(pw)+ 



KN-n)/2] 



k Q2k-l 



k=l 



-(n) 

^2fc dt 2k ~ 1 ' k ' 



if 3 = (-l)V fc T ( 4 l) A k (P (n+2k) ) + ^B^X^M 2 ^ 21 - 1 ^^ 2 ^) ) , (A5) 

1=0 



k-l 



and 



M 



(n) 



[(JV-n-l)/2] Q2fc fc 



t(m^) +r [ ]T ^_^c(;) A ^ 2fc ^(M( n+2fe )) 



fc=i 



The coefficients are given in a compact algebraic form in Ref. 



9. 



(A6) 



4"> 
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B (n) 

c 



2 k c 2k n + 2k'' 
r_l\k-i+i. 



at 



n(n + 2l)\ 



2 i c 2k+2 ( n + 2k + l)(n + 2k + l)\ 
(») . . (-I)*"' n(n + 2l)\ 



k,l 



(A7) 



2 l c 2k {n + 2k)(n + 2k)\' 

The reader is encouraged to apply these formulae for the cases N > 4 and to find convenient 
intermediary calculation. We point out that in equation (lA5h the normalization of the 
quantities is chosen such that the quantities coincide with the electric toroidal 
moments given in literature at least for n up to 3. We present the results given in Ref.- 
concerning the total radiated power X, the recoil force F r and angular momentum loss 
dL/dt: 



1 



Q- 



n + 1 



47re c 3 ^ nn\(2n + l)!!c 2n 

n>l 



(M >n+ i.M in+1 ) + -jP^.P^ 



a' 



(A8) 



/'() 



y- 

/ j r 2n 



U + 2 r~(n) ~ C 2 ,~(n) p(n+l)-, 

:3^c 2 - t(n + l)!(2n + 3)!! L n+1 ' n+2 + a 2 1 " n+2 J 



(n) 



a n!n(2n + 1)!! 
with the obvious notation, 



e i e„ fe (M n 



+i 



.??. 



~(n) . 



(A9) 



■i n - 1 jBi 1 ...i n _ 1 k 
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P „ -n 



n 



an\(2n + 3)\\ 



1 P 

r(n+l) 
r ,n+l 



M 



(n) 
,n+l 



l.M 



(n) 
,n+l 



M 



(n+l) ~(n) 



,n+l 



,n+l 
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in) 



+ 



.("+1) 
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